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Abstract 

We present a new class of asymptotically flat charge static solutions in third order Lovelock 
gravity. These solutions present black hole solutions with two inner and outer event horizons, 
extreme black holes or naked singularities provided the parameters of the solutions are chosen 
suitable. We find that the uncharged asymptotically flat solutions can present black hole with two 
inner and outer horizons. This kind of solution does not exist in Einstein or Gauss-Bonnet gravity, 
and it is a special effect in third order Lovelock gravity. We compute temperature, entropy, charge, 
electric potential and mass of the black hole solutions, and find that these quantities satisfy the 
first law of thermodynamics. We also perform a stability analysis by computing the determinant 
of Hessian matrix of the mass with respect to its thermodynamic variables in both the canonical 
and the grand-canonical ensembles, and show that there exists only an intermediate stable phase. 
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I. INTRODUCTION 



Both string theory as well as brane world cosmology assume that spacetime possesses 
more than four dimensions. In string theory, extra dimensions were promoted from an inter- 
esting curiosity to a theoretical necessity since superstring theory requires a ten-dimensional 
spacetime to be consistent from the quantum point of view, while in the brane world cos- 
mology (which is also consistent with string theory) the matter and gauge interaction are 
localized on a 3-brane, embedded into a higher dimensional spacetime and gravity propa- 
gates in the whole of spacetime. This underscores the need to consider gravity in higher 
dimensions. In this context one may use another consistent theory of gravity in any di- 
mension with a more general action. This action may be written, for example, through the 
use of string theory. The effect of string theory on classical gravitational physics is usually 
investigated by means of a low energy effective action which describes gravity at the classical 
level This effective action consists of the Einstein-Hilbert action plus curvature-squared 
terms and higher powers as well, and in general give rise to fourth order field equations and 
bring in ghosts. However, if the effective action contains the higher powers of curvature 
in particular combinations, then only second order field equations are produced and conse- 
quently no ghosts arise Q. The effective action obtained by this argument is precisely of 
the form proposed by Lovelock j^: 

, m\ 

I G = d d x^/^j^2a k £ k (1) 

J k=0 

where [z] denotes integer part of z, a k is an arbitrary constant and C k is the Euler density 
of a 2fc-dimensional manifold, 

In Eq. (J2J) 8pla\...pua k ls ^ ne generalized totally anti-symmetric Kronecker delta and R^J** * s 
the Riemann tensor. It is worthwhile to mention that in d dimensions, all terms for which 
k > [d/2] are identically equal to zero, and the term k = d/2 is a topological term. So, only 
terms for which k < d/2 are contributing to the field equations. 

In this paper we want to restrict ourself to the first four terms of Lovelock gravity. The 
first term is the cosmological term which we ignore it in the investigation of the proper- 
ties and thermodynamics of the solutions. The second term is the Einstein term, and the 



third and fourth terms are the second order Lovelock (Gauss-Bonnet) and third order Love- 
lock terms respectively. From a geometric point of view, the combination of these terms in 
seven- dimensional spacetimes, is the most general Lagrangian producing second order field 
equations, as in the four- dimensional gravity which the Einstein-Hilbert action is the most 
general Lagrangian producing second order field equations. Here, we will obtain asymptot- 
ically flat solutions of third order Lovelock gravity and investigate their thermodynamics. 

Indeed, it is interesting to explore black holes in generalized gravity theories in order to 
discover which properties are peculiar to Einstein's gravity, and which are robust features of 
all generally covariant theories of gravity. Due to the nonlinearity of the field equations, it is 
very difficult to find out nontrivial exact analytical solutions of Einstein's equation with the 
higher curvature terms. In most cases, one has to adopt some approximation methods or 
find solutions numerically. However, exact static spherical 



of the Gauss-Bonnet gravity have been found in Ref. "\ _ 
Gauss-Bonnet and Einstein-Born-Infeld-Gauss-Bonnet models in Refs. 
solutions with nontrivial topology in this theory have been also studied in Refs. [£ 



symmetric black hole solutions 
and of the Einstein-Maxwell- 
ack hole 
Ej. The 

thermodynamics of the uncharged static spherically black hole solutions has been considered 
in and of charged solutions in [?], [nj]. All of these known solutions in Gauss-Bonnet 
gravity are static. Recently one of us has introduced two new classes of rotating solutions 
of second order Lovelock gravity and investigate their thermodynamics 



Also the static spherically symmetric solutions of the dimensionality continued gravity 
have been explored in Ref. [13j, while black hole solutions with nontrivial topology in this 
theory have been studied in Ref . "\\ . The thermodynamics of these solutions have been 



investigated in Refs. 
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1^ | . Up to our knowledge, no asymptotically flat solution for 
Lovelock gravity higher than second order (Gauss-Bonnet) has been obtained till now with 
two or more fundamental constants. Indeed, the asymptotically flat solution of Ref. |l6j has 
only one fundamental constant which is the gravitational constant G. In this paper we want 
to find new static solutions of third order Lovelock gravity which are asymptotically flat and 
contain two and three fundamental constants and investigate their thermodynamics. As we 
will show later, these asymptotically flat solutions have some properties which do not occur 
in Einstein or Gauss-Bonnet gravity. 

The outline of our paper is as follows. We give a brief review of the field equations of 
third order Lovelock gravity in Sec. |Hj In Sec. we present the static solutions of third 
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order Lovelock gravity in the presence of electromagnetic field with special values of a 2 
and a3, and investigate their properties. In Sec. |lV|we obtain mass, entropy, temperature, 
charge, and electric potential of the d-dimensional black hole solutions and show that these 
quantities satisfy the first law of thermodynamics. We also perform a local stability analysis 
of the black holes in the canonical and grand canonical ensembles. In Sec. |VJ we introduce 
the general asymptotically flat solutions with three fundamental constants and investigate 
their thermodynamics. We finish our paper with some concluding remarks. 

II. FIELD EQUATIONS 

The main fundamental assumptions in standard general relativity are the requirements 
of general covariance and that the field equations for the metric be second order. Based 
on the same principles, the Lovelock Lagrangian is the most general Lagrangian in classical 
gravity which produces second order field equations for the metric. The action of third order 
Lovelock gravity in the presence of electromagnetic field may be written as 



where A is the cosmo logical constant, a 2 and a 3 are Gauss-Bonnet and third order lovelock 
coefficients, C\ = R is just the Einstein-Hilbert Lagrangian, £ 2 = R^wysR^ —^R^uR 111 '+R 2 
is the Gauss-Bonnet Lagrangian, and 



f o Tyfiua-K td opr , o TDP U T} aK DP T I 0/1 DM 1 " 7 " E> DP 

l 3 — ztr tx aR p T tx M „ + up ti VT tir pK + i^ar ti aKvp nr ^ 

+3RR^ K R aKllu + 2AR^ K R ap R KU + 1QR^R U(7 R\ - URR^R^ + R 3 (4) 



is the third order Lovelock Lagrangian. In Eq. Q = d fl A u — d u A p is electromagnetic 
tensor field and is the vector potential. 

Since in Lovelock gravity, only terms for which k < d/2 are contributing to the field 
equations and we want to consider the third order lovelock gravity, therefore we consider the 
(i-dimensional spacetimes with d > 7. Varying the action with respect to the metric tensor 
g^v and electromagnetic tensor field the equations of gravitation and electromagnetic 
fields are obtained as: 





(5) 
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V V F" V = (6) 

where T pu = 2F P IM F pi/ — \F pa F pa g pv is the energy-momentum tensor of electromagnetic field, 
G$ is just the Einstein tensor, and G$ and G$ are given as 

G$ = 2(— R p(TKT R KTa v — 2R ppva R pa — 2R pa R a v + RR pu ) — -£29^ 

/^(3) Of A -QTpCFK p T>\ OD r P P CTK P^ I 9 P TffK f? P^ 

^/ii/ ~~ o^n n aK \ p n VTfl on Acr n T/i n ypK -r ^-n-j, n aK \ p rL T ^ 

T>TpoK p P I Q DT DCK pp I Q DC DTO p/c 

RanrpRup + S>R V(jp R rfj, K + ° H vtk R ap R p 

I A p ttk p pp A p T<7K p pp _|_ T}TpCFK p p I O P P KT P P 

-T t ±ix v ix aKpp rh T rt aKrp rh ~t~ ijl n aKTp ix l/p -|- £ixix v n TpKp 

+% RT vpp RP a RU r ~ % Ra vrp RT a R p ~ ^a^T^P ~ ARR ' ' vpp R ^ \ 

+AR Tp R pT R vp — 8R T u R rp R p p + ARRvpR 11 ^ — R 2 R up ) — -Cj,g pv 



III. STATIC SOLUTIONS 



The metric of <i-dimensional static spherically symmetric spacetime and the vector po- 
tential may be written as: 

ds 2 = -f( r )dt 2 + ^r- + r 2 dn 2 

fin 

K = h(r)6l (7) 

where dfl 2 is the metric of a (d — 2)-sphere. The functions f(r) and h(r) may be obtained 
by solving the field equations (0) and ©. 



A. Seven-dimensional Solutions 



As stated before, the gravitational field equation of third order Lovelock gravity in seven 
dimensions is the most general second order differential equation which presents the solu- 
tions of gravity. Indeed, the solution of third order Lovelock gravity in seven dimensions is 
the most general solution of gravity, based on the principle of standard general relativity. 
Therefore, first we obtain the seven- dimensional solutions of third order Lovelock gravity 
in the presence of electromagnetic field and investigate their properties. Using Eq. (jHJ) one 
can show that h(r) = g/(4r 4 ), where q is an arbitrary real constant which is related to the 
charge of the solution. To find the function f(r), one may use any components of Eq. 



5 



The simplest equation is the rr-component of these equations which can be written as 

-3os(/ - l) 2 + a 2 (f - l)r 2 - ^ r 5 f + a 2 r\f - if - £(/ - 1) - A r io = |! (8) 

where prime denotes the derivative with respect to r. We consider the solutions of Eq. 
(JHJ) for any arbitrary values of a 2 and «3 in Sec. Here, we study the special case of 
a?3 = 2«2 2 = a 2 / 72. Even for this special case, there exist three fundamental constants, A, 
G = 1 and a in the solution, while the solution of ref. Q| has two fundamental constants. 
As we will see below this solution has some properties which will not happen in the Gauss- 
Bonnet gravity with three fundamental constants. The real solution of Eq. (JHJ) with the 
above values of a 2 and is 

„, s r 2 L A Act 3am 3ag ,2 \ 1 ^ 3 | 

/(r) = 1 + _Jl_^l + _ + ___£.j | (9) 

In the above equations m is an integration constant which is related to the mass of the 
solution. Unlike the solutions in Gauss-Bonnet gravity which have two branches, here the 
solution (JHJ) has only one branch. Indeed, Eq. (JHJ) with the above a's has the real solution 
(jnj) and two complex solutions which are the complex conjugate of each other. This feature 
is the same as the asymptotically anti-de Sitter (AdS) solution of Ref. which has a 
unique anti de Sitter vacuum. 

In order to study the general structure of this solution, we first look for the asymptotic 
behavior of the solutions. It is easy to find out that these solutions are asymptotically flat 
for A = 0, AdS for A < and de Sitter (dS) for A > 0. In this paper we are interested in 
the case of asymptotically flat solutions, and therefore we put A = 0. One can show that 
the Kretschmann scalar R fJ , l/ \ K R lluXK diverges at r = 0, and therefore there is a curvature 
singularity located at r = 0. Now we look for the existence of horizons, and therefore we 
look for possible black hole solutions. The horizons, if any exist, are given by the zeros of 
the function f(r) = g rr . 

First, we consider the uncharged solutions, where the horizon(s) is located at 

r+ = { ~ o + ^ ( 10 ) 



2 30 

\j 2 

a VlSOm — 75a 2 1 

2 30 J ^ ' 

As one can see from Eqs. (jlOJ) . and (|llj). for a > there exists only one horizon provided 
the mass parameter, m, is greater than m cr it = 5/3a 2 . That is, there exist a minimum 



value for the mass (or radius of horizon) in order to have uncharged asymptotically flat 
black hole solution. More interesting is the uncharged solutions of third order Lovelock 
gravity with a < 0. In this case there exist an extreme value for the mass parameter, 
m ext = 5/12a 2 . The uncharged solution presents a black hole with inner and outer horizon 
provided m > m ext , an extreme black hole for m = m ext and a naked singularity otherwise. 
This property happens only for third order Lovelock gravity and does not occur in Einstein 
or Gauss-Bonnet gravity. Thus, it is natural to suppose that new solutions in higher order 
Lovelock gravity might provide us with a new window on some corners of low energy limit 
of string theory. 

Second, we consider the charged solutions. It is easy to show that the solution presents 
a black hole solution with two inner and outer horizon, provided the charge parameter q is 
less than, g ext , an extreme black hole for q = q ext and a naked singularity otherwise, where 
<?ext is 

q ext = ^ 1 5760m 2 + 13200a 2 m - 7625a 4 + 5av/l5(96m - 25a 2 ) 3 } (12) 



B. (i-dimensional Solutions 



The d- dimensional static solutions in third order Lovelock gravity may be obtained by 
solving Eqs. (jSj) and (JSJ for the metric given in Eq. (|7jl. Using Eq. one can show that 
h(r) = q/[(d — 3)r d ~ 3 }, where q is an arbitrary real constant which is related to the charge 
of the solution, and the function f(r) is the solution of the following equation: 

d-2 



180(V)«ar(/ - I) 2 - 6(V)a 2 (/ - l ) r + 



f + Ar 6 



5 J»*>KJ V "V 3 ™ ■ 2 

+360(t 2 )a 3 (/ - l) 3 - 12( d 4 2 )a 2 r 2 (/ - l) 2 + ( d " 2 )r 4 (/ - 1) = -gV°- 2d (13) 

where prime denotes the derivative with respect to r. The solution of Eq. (fT3*|) for arbitrary 
values of a 2 and a 3 will be introduced in Sec. El Here we consider the solutions of Eq. (TH?j) 
for aj's given as 



a 



a 



a 3 



a 2 



72( d " 3 )' ^ (d-3)(d-4) 



(14) 



Again, Eq. (|T3*j) with condition ()14|) has one real and two complex solutions which are the 
complex conjugate of each other. The real solution of Eq. (fT3*j) with condition (fT4*j) is 

2 ' / 6Aa 3am 6ag 2 x ' ' 



f(r) 



1 + 



1 



a 



1 + (d - 1) (d - 2) + (d - 2)r rf "! (d - 2) (d - 3)r M " 4 



(15) 



where m is the mass parameter. 

These solutions are asymptotically flat for A = and AdS or dS for negative or positive 
values of A respectively. Also one can show that the Kretschmann scalar diverges at r = 
and therefore there is a curvature singularity located at r = 0. 

In order to investigate the existence of black hole solutions for A = we, first, consider 
the uncharged solutions. As in the case of seven-dimensional solutions, one can show that 
for a > 0, there exists a critical value of horizon radius, r +cr j t (or a critical value of mass, 
m crit ) such that the solution presents a black hole for r + > r +crit , and a naked singularity 
otherwise, where r +crit is the larger real root of the following equation: 

3(d - 3)rl + 3{d - 5)ar 2 + + (d- 7)a 2 = (16) 

The critical mass, m CTit now may be obtained easily as 

(d - 2) \d + 3 + y/-3(d- l)(d- 9) 



m 



crit 



3d - 15 - y/-3(d- l)(d-9) 



( d - 7 )/ 2 , \ (d-3)/2 



a 



y{ ' ] {—s) < 17 > 

Note that the critical mass is real only for d < 9, and therefore there exist no uncharged 
black hole in dimension greater than nine for «2 and 0:3 given in Eqs. (j!4jl . Unlike to this 
special case, we will see in Sec. El that black hole solutions exist in any dimension for the 
arbitrary values of «2 and a^. 

For a < 0, the ci-dimensional solutions of third order Lovelock gravity present a new 
property that does not occur in lower order of Lovelock gravity. Indeed the uncharged 
solution for negative values of a presents a black hole with two inner and outer horizon, 
provided r + > r +ext , an extreme black hole for r + = r +ex t, and a naked singularity otherwise, 
where r +ext is the smaller solution of Eq. 1)1 7]). In this case the extremal value of mass is 



(d - 2) d + 3 - y/-3{d- l)(d-9) 

fT^ext. = Z ' 



3d-lh+ ^-3{d- l)(d-9) 



{d - 1)/2 / I I \ (d"3)/2 



6 Vrf-3 



(18) 



Second, we consider the charged solutions. It is easy to show that the solution presents a 
black hole solution with two inner and outer horizon, provided q is less than q ext , an extreme 



8 



black hole for q = q ext , and a naked singularity otherwise , where q ext is the solution of 

3(d - 5)ar 2 + d - 8 + 3(d - 3)r^" 6 + (d - 7)a 2 r 2 + d - 10 - Q(d - 2)"Y = (19) 

IV. THERMODYNAMICS OF BLACK HOLES 

One can obtain the temperature of the event horizon by analytic continuation of the 
metric. The analytical continuation of the Lorentzian metric by t —>■ ir yields the Euclidean 
section, whose regularity at r = r + requires that we should identify r ~ r + f3 + , where (3+ 
is the inverse Hawking temperature of the horizon given as 

x f'(r) _ (d - 2)rl d - w [3(d - 3)rj + 3(d - 5)ar 2 + + (d - 7)a 2 } - 2q 2 

/+ + 4vr 127r(r* + a)V^ 9 1 J 

Usually entropy of the black holes satisfies the so-called area law of entropy which states 
that the black hole entropy equals to one-quarter of horizon area One of the surprising 
and impressive feature of this area law of entropy is its universality. It applies to all kind of 
black holes and black strings of Einstein gravity jscj . However, in higher derivative gravity 



the area law of entropy is not satisfied in general 
Lovelock gravity is [22] 

( [(d-l)/2] 

s 



2l\ . It is known that the entropy in 



- ka * I d d - 2 x^C k ^ (21) 

fc=i J 

where the integration is done on the (d — 2)-dimensional spacelike hypersurface of Killing 
horizon, g^ v is the induced metric on it, g is the determinant of g^ v and is the fcth order 
Lovelock Lagrangian of g^ u . Thus, the entropy in third order Lovelock gravity is 

S = l - J d d - 2 xyfg~ (l + 2a 2 R + 3a z (R^ UK R^ K - R^BT + i? 2 )) (22) 

where R^ vpa and R pu are Riemann and Ricci tensors and R is the Ricci scalar for the induced 
metric g a b on the (d — 2)-dimensional horizon. It is a matter of calculation to show that the 
entropy of black holes is 



S^ i J 1 + ^ + ^ ] (23) 



where S^_2 is the area of a unit radius {d — 2)-dimensional sphere. 

The charge of the black hole can be found by calculating the flux of the electric field at 
infinity, yielding 



The electric potential $, measured at infinity with respect to the horizon, is defined by 
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$ = A mX m Uoo-^X M U + > (25) 
where x = d/dt is the null generator of the horizon. One finds 

$ = 7 Wl (26) 

(d - 3)r d + - 3 V ; 

The mass of black hole can be obtained by using the behavior of the metric at large r. 
It is easy to show that the mass of black hole is 

(d -2)(d- 3)(3rj + Sarj + c?)r^ + 2g 2 
16n m ~ (rf-3)rt 3 10 

We now investigate the first law of thermodynamics. Using the expression for the entropy, 
the charge and the mass given in Eqs. (|23|). (|24|) and ()27j) . one can compute dM/dr, dS/dr 
and dQ/dr. Then, by using the chain rule, it is easy to show that the quantities 

T =( — ) , $= l^» 5 (28) 



dSJ' \dQ 



are exactly the same as the temperature and electric potential given in Eqs. (J2~U|) and (f2*U|) 
respectively. Thus, the thermodynamic quantities calculated in Eqs. (j2*Uj) and satisfy 
the first law of thermodynamics, 

rfM = TciS + <$>dQ (29) 



A. Stability in the canonical and the grand-canonical ensemble 

The stability of a thermodynamic system with respect to the small variations of the 
thermodynamic coordinates, is usually performed by analyzing the behavior of the entropy 
S(M, Q) around the equilibrium. The local stability in any ensemble requires that S(M, Q) 
be a convex function of their extensive variables or its Legendre transformation must be 
a concave function of their intensive variables. Thus, the local stability can in principle 
be carried out by finding the determinant of the Hessian matrix of S with respect to its 
extensive variables, Hjr iJf - = [d 2 S/dXidXj], where X^s are the thermodynamic variables 
of the system. Indeed, the system is locally stable if the determinant of Hessian matrix 

n 

satisfies x , < |2J. Also, one can perform the stability analysis through the use of the 
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determinant of Hessian matrix of the energy with respect to its thermodynamic variables, 
and the stability requirement H^. x . < may be rephrased as H^y. > 25]. 

The number of the thermodynamic variables depends on the ensemble which is used. In 
the canonical ensemble, the charge is a fixed parameter, and therefore the positivity of the 
heat capacity Cq = T(dS/ dT)q is sufficient to assure the local stability. The heat capacity 
for the solutions given by Eq. ()15|) shows that there exist an upper limit for the radius of 
horizon, r +st which is the real root of St can = 0, where St can is 

Stcan = (d - 2)[(d - 7)a 3 + 2r +st \d - 10)a 2 - 18r +st 4 a + 3r +st 6 (rf - 3)] 

-6(2d - 5)r +s r M+12 g 2 - 6(2d - 9)r +st - 2d+10 aq 2 (30) 

Therefore, the charged solutions are stable provided the horizon radius lies in the range 
r+ext < r + < r +st, where r +ext is the solution of Eq. (fT9j) . and r +st is the real solution of Eq. 
(|3()j) . Indeed, for r + < r +cxt , we have no black hole while for r + > r +st , the black hole is not 
stable and therefore we have only an intermediate stable phase. Note that for the case of 
uncharged black holes and a > 0, r + should be greater than r +crit as stated before. Thus 
for positive a, the horizon radius of stable black holes lies in the range r +cr i t < r + < r +st , 
where r +crit is the larger root of Eq. 1)1 7]). Of course, for negative a, which one can have 
a black hole with inner and outer horizon, the black hole is stable provided r + lies in the 
range r +ext < r + < r +st , where now r +ext is the smaller root of Eq. (fTTj) , and again we have 
only an intermediate stable phase. 

In the grand-canonical ensemble, the stability analysis can be carried out by calculating 
the determinant of Hessian matrix of the energy with respect to S and Q. The zeros of the 
determinant of Hessian matrix are given by St gc = 0, where St gc is 

St gc = (d - 2)[(d - 7)a 3 + 2 r 2 +st (d -10) a 2 -18 r 4 +st a + 3 r 6 +st (d- 3)] 

+18r;^V«-6r; s 2 t ^V (31) 

Three cases happen for the roots of St gc = 0: 

1. It has two real roots r +st i < r +st 2. In this case, the determinant of Hessian matrix is 
positive provided the radius of horizon lies in the range between r + i ar g er < r + < r +st 2, where 
f+iarger is the larger value between r +ext and r +st i. Thus, there exist an intermediate stable 
phase 

2. It has one real root r +st2 . In this case, there exists an intermediate stable phase with 
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radius of horizon between r +cxt and r +st2 . 

3. It has no real solution. In this case the black hole is unstable for the whole range of r+. 

Numerical analysis shows that r +st2 < r +st (if r +st2 exists), and therefore the region of 
stability is smaller for the grand-canonical ensemble. This is due to the fact that the number 
of thermodynamic variables in the canonical ensemble is less than that of the grand-canonical 
ensemble. 



V. GENERAL SOLUTIONS 



Finally we give the general solutions of third order Lovelock gravity in d dimensions for 
any arbitrary values of ot 2 and a 3 . In this case the solution of Eq. (fT^j) is 

2(d+l) /in 1 \ 

,(f) ( a, ^1/3 r 3 (73=^ -57^*3) 
/(r) = l + r 2 ^ + - (b+^y/b) + — ^TtV^ (32) 



C2 5 ) «3 V- 24 *7 36a 3 (^ + ffv^) 1/3 

where ^ and Q are 

a^r 2 ^ 2 ) a 2 a 3 r 2( - d ^ ma\(d — 2)r d ~ 3 a|g 2 



216( d 2 5 ) 3 144( d - 3 )(^-6)(rf-5) 120( d - 2 ) 120( d ~ 2 )(d - 3) 



a 3 r 4 ^ a 2 r 4 ^ 2 1± ^a 2 mq 2 (d - 2)r d ' 3 

648( d J 3 ) 3 ~ 36(rf - 6) 2 (d - 5) 2 ( d J 3 ) 2 + 25(d - 3)( d ~ 2 ) 2 
a 2 a 3 m{d - 2)r 3d ~ 7 ^ a 2 m 2 (d - 2) 2 r 2 ( d - 3 ) 



15 x 2 1 -(- 1 ) d (rf - 6)(d - 5)( d 5 2 )( d J 3 ) 25 x 4M-i) d ( d s 2 ) 2 



a 2 a 3 q 2 r 2 ^ 8 x 2 i±i ^ (d - 2)ma 3 r M " 7 



+ 



3(t 3 )( d 4 2 )( rf - 3)(d- 5)(d - 6) 2 45(d - 5) 3 (d - 6) 3 ( d ~ 2 ) 
q 4 a\ 16g 2 a 3 r 2 ( d - 2 ) 



25(d - 3) 2 (^ 2 ) 2 9(d - 3)(d - 2)(d - 6) 3 (d - 5) 3 (V) 
In the above equations m is an integration constant which is related to the mass of the 
solution. This solution is asymptotically flat and has a curvature singularity at r = 0. 
Note that this asymptotically flat solution has three fundamental constants G = 1, a 2 and 
0:3. Numerical analysis shows that the solutions (|32J) may present black holes with two 
inner and outer horizons, extreme black holes, black holes with one event horizon, or naked 
singularities depending on the values of a 2 , a 3 , m and q. To be more clear we give the range 
of values of m, a 2 and a 3 for the 7-dimensional uncharged solution in order to have black 
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FIG. 1: m versus «2 for q = 0, d = 7, and 03 = 3«2/2. 



1- 




FIG. 2: f(r) versus r for d = 10, q = 0, m = 500, 02 = 1-0, and 03 = 2.0. 

hole solutions. For 0:3 = 3a2/2, any solution with m and 0:2 in region I of Fig. ^represents 
a black hole solution with an event horizon. 

Unlike to the case of solutions with special values of and given in Sec. 11111 the 
uncharged solution for arbitrary values of fundamental constants «2 and a.3 can have horizon 
even for d > 9. Figure El shows the function f(r) for d = 10 as a function of r. It shows that 
there is an event horizon at r = 2. 

To investigate the first law of thermodynamics for the solutions with arbitrary values of 
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«2 and «3, one should perform a numerical analysis. Using the expression for the entropy 
and the charge given in Eqs. (|2*3j) and ([24)1 . and the fact that M = (167r)]m, where m 

is the solution of f(m,r + ) = 0, one can show numerically that the temperature and electric 
potential of Eq. (J28)) are exactly equal to the temperature and electric potential which can 
be computed through the use of geometrical analysis. Thus, the first law of thermodynamics 
is valid for the general solution given in Eq. (JH2*j) . Also, one can perform a stability analysis 
in canonical and grand-canonical ensembles numerically. Again, as in the case of special 
solutions of third order Lovelock gravity, numerical analysis shows that there exist only 
an intermediate stable phase for black hole solutions which is smaller in grand-canonical 
ensemble . 

VI. CLOSING REMARKS 

In this paper, we added the third order Lovelock terms to the Gauss-Bonnet-Maxwell 
action, and introduced a new class of static solutions which are asymptotically flat, AdS or 
dS for A = 0, A < or A > respectively. We were only interested in the asymptotically 
flat solutions. 

First, we consider the solutions for special values of and 03 given in Eq. (|14|). For 
uncharged solutions, we found that for a > 0, the solutions present black holes provided 
the mass is larger than a critical value, m cr it, given in Eq. (jlfij) . For negative a, they 
present black holes with two inner and outer horizons for m > m ext , extreme black holes for 
m = m ext or naked singularity for m < m ext , where m ext is given in Eq. (jl8j) . These kind of 
uncharged solutions exist only in third order Lovelock gravity and do not happen in Einstein 
or Gauss-Bonnet gravity. For the case of charged solutions, we found that there exists an 
extremal value of r + , given by Eq. (|19|). that determines whether the solution presents a 
black holes with two inner and outer horizons, an extreme black hole or a naked singularity. 
Accordingly we obtained temperature, entropy, charge , electric potential and mass of these 
black hole solutions. We also investigated the first law of thermodynamics, and found 
that these thermodynamics quantities satisfy the first law of thermodynamics. Also, we 
performed a stability analysis in canonical ensemble by considering the heat capacity of the 
solution and found that there exist two unstable phases separated by an intermediate stable 
phase. This analysis was also done through the use of the determinant of the Hessian matrix 
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of M(S, Q) with respect to its extensive variables and we got the same phase behavior with 
a smaller region of stability in the grand-canonical ensemble or no stable phase depending 
on the values of black hole's parameters. 

Second, we introduced the general solution of third order Lovelock gravity for arbitrary 
values of «2 and 0:3 and investigate its properties. Numerical analysis showed that these 
solutions may be interpreted as black hole solutions with two inner and outer event horizons, 
extreme black holes or naked singularity depending on the parameters of the solutions. We 
also found that the conserved and thermodynamics quantities for these general solutions 
satisfy first law of thermodynamics. Again, the stability analysis showed that there exists 
only an intermediate stable phase, for the black hole solutions which is smaller for the 
grand-canonical ensemble. 

As we motion, the as ym pto tl ca„ y flat 8 o. Uti o„ 8 obtained in Q contain onl y one funda- 
mental constant. Finding new solutions in continued Lovelock gravity with more fundamen- 
tal constants remains to be carried out in future. Also, the generalization of these solutions 
to the case of rotating solutions will be given elsewhere. 



[1] M. B. Greens, J. H. Schwarz and E. Witten, Superstring Theory, (Cambridge University Press, 
Cambridge, England, 1987); D. Lust and S. Theusen, Lectures on String Theory, (Springer, 
Berlin, 1989); J. Polchinski, String Theory, (Cambridge University Press, Cambridge, England, 
1998). 

[2] B. Zwiebach, Phys. Lett. B156, 315 (1985); B. Zumino, Phys. Rep. 137, 109 (1986). 

[3] D. Lovelock, J. Math. Phys. 12, 498 (1971); N. Deruelle and L. Farina-Busto, Phys. Rev. D 

41, 3696 (1990); G. A. MenaMarugan, ibid. 46, 4320 (1992); 4340 (1992). 
[4] D. G. Boulware and S. Deser, Phys. Rev. Lett., 55, 2656 (1985); J. T. Wheeler, Nucl. Phys. 

B268, 737 (1986). 
[5] J. T. Wheeler, Nucl. Phys. B273, 732 (1986). 
[6] D. L. Wiltshire, Phys. Lett. B 169, 36 (1986). 
[7] D. L. Wiltshire, Phys. Rev. D 38, 2445 (1988). 

[8] R. G. Cai and K. S. Soh, Phys. Rev. D 59, 044013 (1999); R. G. Cai, ibid. 65, 084014 (2002). 
[9] Y. M. Cho and I. P. Neupane, Phys. Rev. D 66, 024044 (2002). 



15 



[10] R. C. Myers and J. S. Simon, Phys. Rev. D 38, 2434 (1988); R. C. Myers, Nucl. Phys. B289, 
701 (1987). 

[11] J. E. Lidsey, S. Nojiri and S. D. Odintsov, J. High Energy Phys., 06, 026 (2002); M. Cvetic, 
S. Nojiri, S. D. Odintsov, Nucl. Phys. B628, 295 (2002); S. Nojiri and S. D. Odintsov, Phys. 
Lett. B 521, 87 (2001); ibid. 523, 165 (2001). 

[12] M. H. Dehghani, Phys. Rev. D 67, 064017 (2003); ibid. 69, 064024 (2004); ibid. 70, 064019 
(2004). 

[13] M. Banados, C. Teitelboim and J. Zanelli, Phys. Rev. D 49, 975 (1994). 

[14] R. Aros, R. Troncoso and J. Zanelli, Phys. Rev. D 63, 084015 (2001). 

[15] J. P. Muniain and D. Piriz, Phys. Rev. D 53, 816 (1996). 

[16] J. Crisostomo, R. Troncoso and J. Zanelli, Phys. Rev. D 62, 084013 (2000). 

[17] T. Clunan, S. F. Ross and D. J. Smith, Class. Quant. Grav. 21, 3447 (2004). 

[18] F. Muller-Hoissen, Phys. Lett. B 163, 106 (1985). 

[19] J. D. Beckenstein, Phys. Rev. D 7, 2333 (1973); S. W. Hawking, Nature (London) 248, 30 

(1974); G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2738 (1977). 
[20] C. J. Hunter, Phys. Rev. D 59, 024009 (1999); S. W. Hawking, C. J. Hunter and D. N. Page, 

ibid. 59, 044033 (1999); R. B. Mann, ibid. 60, 104047 (1999); ibid. 61, 084013 (2000); 
[21] M. Lu and M. B. Wise, Phys. Rev. D 47, R3095,(1993); M. Visser, ibid. 48, 583 (1993). 
[22] T. Jacobson and R. C. Myers, Phys. Rev. Lett. 70, 3684 (1993); R. M. Wald, Phys. Rev. D 

48, R3427, (1993); M. Visser, ibid. 48, 5697 (1993); T. Jacobson, G. Kang and R. C. Myers, 

ibid. 49, 6587,(1994); V. Iyer and R. M. Wald, ibid. 50, 846 (1994). 
[23] M. M. Caldarelli, G. Cognola and D. Klemm, Class. Quantum Grav. 17, 399 (2000). 
[24] M. Cvetic and S. S. Gubser, J. High Energy Phys. 04, 024 (1999). 
[25] S. S. Gubser and I. Mitra, J. High Energy Phys., 08, 018 (2001). 



16 



